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Abstract 



The anomalous pseudoscalar-photon vertex is studied in real time in and out 
of equilibrium in a constituent quark model. The goal is to understand the 
in-medium modifications of this vertex, exploring the possibility of enhanced 
isospin breaking by electromagnetic effects as well as the formation of neutral 
pion condensates in a rapid chiral phase transition in peripheral, ultrarela- 
tivistic heavy-ion collisions. In equilibrium the effective vertex is afflicted by 
infrared and pinch singularities that require hard thermal loop (HTL) and 
width corrections of the quark propagator. The resummed effective equilib- 
rium vertex vanishes near the chiral transition in the chiral limit. In a strongly 
out of equilibrium chiral phase transition we find that the chiral condensate 
drastically modifies the quark propagators and the effective vertex. The ensu- 
ing dynamics for the neutral pion results in a potential enhancement of isospin 
breaking and the formation of 7r° condensates. While the anomaly equation 
and the axial Ward identity are not modified by the medium in or out of 
equilibrium, the effective real-time pseudoscalar-photon vertex is sensitive to 
low energy physics. 
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I. INTRODUCTION 



The forthcoming ultrarelativistic heavy ion colliders, RHIC at BNL and LHC at CERN 
will allow us to probe the physics of the quark gluon plasma, the hadronization and the chiral 
phase transitions, providing for the first time, experimental access to a phase transition in 
particle physics as predicted by QCD Current experiments at CERN SPS are seeking 
manifestations of these phase transitions with the large acceptance detectors NA49 and 
WA98, while the large acceptance detector capabilities at RHIC will allow an event-by- 
event analysis of both hadronic and electromagnetic signals. Currently WA98 at SPS-CERN 
studies fluctuations in the neutral to charged pion ratios on an event-by-event basis [TJ, 
and previously the MINIMAX collaboration at Fermilab studied similar signals [[3]. These 
experiments are particularly important in the detection of disoriented chiral condensates 
(DCC) which are coherent regions of chirally misaligned pion condensates, proposed and 
discussed by many authors |4j] , originally within the context of CENTAURO events in cosmic 
ray experiments || but more recently as a possible signal for the chiral phase transition if it 
occurs out of equilibrium With the turning on of the above mentioned experiments, 

with their capacity to analyze data from ultrarelativistic heavy-ion collisions on an event- 
by-event basis, a theoretical exploration of various signals of the chiral phase transition is in 
order. 

The axial anomaly |IJ has long been a source of extremely interesting physics, ranging 
from the explanation of the 7r° — > 77 decay to understanding the mass splitting between 
the T) and the 7/ pseudoscalar mesons |IT]| . In this article we focus on the equilibrium and 
nonequilibrium aspects of the anomalous electromagnetic interaction of the neutral pion to 
investigate the possibility of experimentally important signals arising from this interaction 
and its modifications in a medium. In particular, since we expect that the collision region 
will contain a large amount of energy density, a finite temperature study as well as a study 
of the effects of the axial anomaly in out-of-equilibrium situations will prove important for 
our understanding of forthcoming experimental results. 

Detailed studies in equilibrium have found that the anomalous Ward identity is inde- 
pendent of temperature and chemical potential [|l^] and that in chiral perturbation theory, 
valid for T << T c the amplitude for neutral pion decay into two photons is independent of 
temperature. 

However, in a series of recent articles, Pisarski [HJ and Pisarski, Trueman and Tytgat 



14] have provided a detailed analysis of the chiral anomaly and the anomalous vertex at 



finite temperature using the Euclidean formulation and concluded, that while the amplitude 
for the neutral pion decay into two photons at zero temperature is directly proportional 
to the coefficient of the axial anomaly (given by the anomalous Ward identity), at finite 
temperatures this relationship no longer holds. This result was anticipated in where it 
was pointed out that at finite temperature there is a new vector (the four-velocity of the 
heat bath) which leads to a more complicated tensor structure for the triangle diagram that 
yields the anomalous amplitude and invalidates the tensor analysis leading to the Sutherland- 
Veltman theorem HTHl. 
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Pisarski [HJ also pointed out that near the critical temperature for the chiral phase 
transition, the amplitude for 7r° — > 77 vanishes. This conclusion was based on an analysis of 
the corresponding 3-point vertex with spacelike external photons, using the imaginary time 



formalism. Gelis |17| studied the anomalous decay of the neutral pion into two photons at 



finite temperature and pointed out some subtleties of the triangle diagram in the limit of 



soft external momenta, and reconciled certain discrepancies between the results of [|13| , |14 | 



and those of fll8| . A study in the linear sigma model at finite temperature and chemical 



potential in equilibrium was anticipated in 19 with similar results to those found in [18,17 



In this work we apply the real time formulation of quantum field theory to the linear 
sigma model (LSM) coupled to constituent quarks and photons and calculate the effective 
coupling between the ir° and photons both at finite temperature as well as out of equilibrium. 
This allows us to obtain far more information from the system than could be found from 
calculations of suitably defined decay amplitudes (for e.g. |I7|); in particular, we will be able 



to construct a real time equation of motion for the neutral pion condensate, and make contact 
with recent proposals of DCC formation p],B-|9j via anomalous electromagnetic interactions 



We first focus our attention on the equilibrium case in the region near the chiral phase 
transition where the pions and constituent quarks are becoming massless. We find that in 
the limit of soft photon momenta near the phase transition two different types of divergences 
appear: infrared divergences in the chirally restored phase associated with the masslessness 
of the quarks in the triangle diagram, and pinch singularities that arise solely from the finite 
temperature contributions. The infrared singularities near the chiral limit are the source of 
discrepancy between the work in |13| , |T4|| performed in the imaginary time formulation and 



that of JT8| in the real time version. This discrepancy and its resolution has been recently 



discussed in |p!7| |. The pinch singularities have a very different origin. For photons of very 



small virtuality, the quarks that are present in the medium and that can contribute to the 
triangle diagram, can propagate on-shell. If the quarks are treated as bare objects, then these 
on-shell states in the thermal bath will propagate without damping and a pinch singularity 
will appear as a consequence of these undamped on-shell intermediate states. 

Both singularities require that the internal quark propagators be dressed by self-energy 
corrections, giving them a chirally symmetric mass as well as a width of thermal origin. In the 
Yukawa theory, the fermion receives Hard Thermal Loop (HTL) corrections to its self-energy 



17| , |2lH when its momentum « gT where g is the Yukawa coupling and T the temperature. 
As long as the quark mass derived from the tree-level Yukawa interaction, m q » g T, HTL's 
are unimportant as far as the real part of the fermion self-energy is concerned. However in 
the critical region in the strict chiral limit, the constituent quark mass m q vanishes and 
the contribution of hard thermal loops (HTL) becomes important. The pinch singularities 
require a further self-energy resummation that incorporates the on-shell width of the thermal 
excitations; this goes beyond the HTL contribution. 

We will argue that both the hard thermal loop correction to the self-energy as well as 
the width are dominated by scalar and pseudoscalar exchange. Furthermore, we find that 
the quark width in the medium arises from the decay of the scalar into quark- ant iquark 
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pairs as found previously in [^] within a different context and from collisional contributions 
in the HTL limit We compute these corrections to the quark propagators explicitly 

in the equilibrium case, and we find, in agreement with |T3|JT^ , p!7|] , that the pion-photon 
vertex vanishes in the strict chiral limit near the critical temperature. 7r° — > 77 at finite 
temperature is, however, of limited relevance for heavy-ion collisions since the typical lifetime 
of the neutral pion is several orders of magnitude larger than that of the fireball (the latter 
being « 50fm/c). The r] meson, however, has a much shorter lifetime ~ 1.2 KeV and 
finite temperature effects may be of importance for the process r] — > 77 which occurs with a 
branching ratio of 40%. Finite temperature effects could also be of importance in the cooling 
of hot and dense stars |18[ , and in certain models of baryogenesis [^3[ . 

Our interest in the nonequilibrium calculation of the pion-photon anomalous vertex comes 
in part from a recent proposal in [M. It is argued there that in the presence of strong 
electromagnetic fields, such as might be found in peripheral collisions between heavy nuclei 
2^|, the anomalous coupling might enhance DCC formation along the the neutral pion 



direction. To quantify the effect of strong EM fields on DCC formation, note that the 
coupling of the neutral pion to the electromagnetic field is given by ~ q^^qF^F^ . In 
peripheral collisions, the strong electromagnetic fields give rise to a semiclassical contribution 
to F^ V F^ U oc E ■ B oc Z 2 e 2 b/R 6 with b the impact parameter |20|1 . Therefore, for large 
impact parameters (very peripheral collisions b >> 2R) and heavy nuclei the anomalous 
interaction of the neutral pion with the electromagnetic field could induce large neutral pion 
condensates. A further interesting scenario was investigated in M, namely the possibility of 
the anomalous vertex inducing photon production via parametric amplification during the 
chiral phase transition. These investigations studied the phenomenological consequences of 
the anomaly by using the vacuum form for the vertex g nll ~ ot/ f n . 

But our previous comments clearly suggest that the vertex is not protected from correc- 
tions by the medium and the strength of the interaction should be drastically modified by 
the nonequilibrium state, as is the case at finite temperature. Our aim therefore is : i) to 
study the nature of these modifications and provide a detailed construction of the vertex 
in real time, and ii) to provide a quantitative assessment of the potential phenomenological 
consequences of these in-medium corrections. In this paper we deal with the first part of 
the program, i.e. the more formal aspects of setting up a framework for studying the out- 
of-equilibrium interactions and relegate to a forthcoming article, the phenomenological and 
numerical implications of these modifications. 

In the nonequilibrium case we assume that the initial state is some thermal state where 
the chiral order parameter has a small amplitude set by the explicit chiral symmetry break- 
ing term. We then assume that the system undergoes a phenomenological quench and the 
order parameter (< o >) then rolls toward the bottom of its potential. In this case the 
effective quark mass is then time dependent which induces a time dependence in the triangle 
diagram and hence in the effective anomalous electromagnetic coupling. Furthermore as the 
< a > evolves, long-wavelength pion fluctuations are generated that feed in to the triangle 
diagram as well. We obtain the quark propagators during the rolling of the chiral conden- 
sate in a mean field approximation. The effective pseudoscalar-photon interactions are then 
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obtained by explicitly integrating out the quarks and computing the triangle diagram in 
real-time and a local limit is extracted for quasistatic classical electromagnetic fields. Using 
this systematic formulation we then derive the equation of motion for the neutral pion field 
after including both (a) self-consistent mean-field effects and (b) effective nonequilibrium 
anomalous couplings to the quasi-static semiclassical electromagnetic fields which are of rel- 
evance in peripheral heavy-ion collisions. We then argue that unstable long-wavelength pion 
fluctuations can enhance the isospin breaking electromagnetic effects via the nonequilibrium 
anomalous coupling. 

In an Appendix we also show that the anomaly equation and the axial Ward identity are 
not modified by the medium in or out of equilibrium. Thus our study shows that the Ward 
identity does not completely determine the vertex in a medium and more importantly, the 
latter is sensitive to low energy physics and is therefore model-dependent. 

In section II we obtain the effective action in real time in equilibrium, analyze the infrared 
divergences and obtain the HTL corrections to the self-energy. We analyze the different 
processes that lead to a width for the fermionic quasiparticles and obtain a quantitative 
estimate for the local limit of the effective pseudoscalar-photon vertex. [By local limit we 
mean taking the external momenta to zero]. 

Section III is devoted to the nonequilibrium case. We consider a 'quenched' chiral phase 
transition in which the expectation value of the chiral order parameter begins with small 
amplitude and rolls down its potential hill towards the minimum equilibrium state. The 
nonequilibrium effective action for pions and photons is obtained in the mean-field approx- 
imation. Finally we obtain the effective nonequilibrium vertex for the neutral pion in the 
limit of quasistatic classical electromagnetic fields, thus obtaining the effective equation of 
motion for the neutral pion including its anomalous coupling. 

Section IV summarizes our results. Appendix A deals with obtaining the anomaly equa- 
tion in real time both in and out of equilibrium, while Appendix B provides the technical 
details of the solution of the Dirac equation in a time dependent scalar background. 

II. 7T77 VERTEX IN EQUILIBRIUM AT FINITE TEMPERATURE 

The purpose of this section is to provide a detailed analysis within real-time, finite tem- 
perature field theory, of the effective anomalous coupling between pseudoscalars and photons 
in a simple model that incorporates all of the essential physics. We work within a Gell-Mann 
Levy type model containing a single fermion flavour ip charged under U(l) em , with Yukawa 
couplings to a scalar o and a pseudoscalar tt. As it will become clear below the results of our 
analysis generalize straightforwardly to the U(2) ® U(2) constituent quark model with the 
appropriate inclusion of the isospin group structure in the vertices. The Lagrangian density 
for the model under consideration is given by 

£ = - rejf)^ - 2g $(a + - ^F^F^ + C n . (2.1) 

C n determines the dynamics of the 7r and a fields and is a 4 scalar theory given by, 
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The vacuum constituent quark mass ~ 300 MeV ~ 2gf n with the vacuum value of 
f n = 93 MeV leads to g ~ 1.5 > e ~ 0.3. Although the large phenomenological value of the 
Yukawa coupling would preclude a perturbative expansion, we will assume 1 >> g » e so 
as to be able to provide a quantitative assessment of the in-medium effects at least within a 
perturbative expansion. While the assumption g << 1 is necessary if we want to provide a 
quantitative understanding based on perturbation theory, the condition e << g is imposed 
only to simplify the calculations and can be easily relaxed without affecting the qualitative 
validity of the results in what follows. 

A small linear term c a ensures that in equilibrium isospin symmetry is broken along the 
a direction. At zero temperature (a) « /„. and the explicit symmetry breaking term gives 
the pions their mass. The strict chiral limit refers to c = 0. At finite temperatures T ^ 0, 
the expectation value of a acquires temperature dependence and we split the quantum field 
into a c-number expectation value and quantum fluctuations: 

a(x) = a (T) + X {x) ; (a{x)) = a (T), (2.3) 

Here the expectation value of the Heisenberg operators is evaluated in the thermal density 
matrix. The fermion thus gets a temperature dependent tree-level mass, 

m q (T) = 2ga (T). (2.4) 

Our subsequent analysis will be valid in the regime where T » g ctq(T). In particular 
this condition will be satisfied near the critical temperature where in the strict chiral limit 
c = 0, <Jq(T) vanishes as a (0) Jl — |r ; T — > T c r. 



A. The real time effective action: 

Within the simple constituent quark model introduced above, we now want to study the 
effective theory of pions interacting with photons at a finite temperature T ^ 0. Moreover, 
we would like to be able to study the dynamics of such interactions. Static phenomena can 
be studied using the imaginary time formulation of finite temperature field theory. 

The proper way to study dynamical phenomena in quantum field theory is via the real 
time evolution of an initially prepared density matrix. This is implemented in the path inte- 
gral formulation by use of the Closed Time Path (CTP) formalism of Keldysh and Schwinger 



25| , |26|| . This formalism allows us to compute the time evolution of expectation values of 
Heisenberg operators and can be applied both to real time equilibrium calculations as well 
as nonequilibrium ones. In this formulation the action is defined on a contour that runs 
forward and backward in time as befits the time evolution of an initially prepared density 
matrix. Field operators on the forward and backward branches are denoted with + and 
— superscripts respectively and are independent field variables. For a description of this 
formulation within the setting of problems akin to those studied here see [E5H29I]. 



5 



As we have stressed in the introduction our goal is to derive an effective description for the 
pions and photons which really means obtaining the real time effective action by integrating 
out the quarks including the anomalous pion-photon interactions. We remark that, for our 



purposes the decay amplitude for n — > 77 as defined in |TjJ for example, is a more restricted 
quantity than the full real-time effective action. The former is given by a particular cut of 
the self-energy graph (Fig. 2) for the 7r. The couplings in the effective action also contain 
discontinuities arising from the on-shell propagation of intermediate fermions in the thermal 
bath. Cross-couplings between fields living on the forward and backward contours will also 
be induced. All of these terms in the effective action will be important for understanding 
the effective pion-photon interactions in the thermal bath; for example the full real-time 
equation of motion for pion condensates will necessarily receive contributions from all these 
different terms. 

Although the effective action in Euclidean time has been derived previously, we obtain 
here the real time effective action, which we believe is a novel approach. The utility of the 
real time effective action lies in the fact that it allows us to obtain the real-time, retarded 
equations of motion for expectation values even out of equilibrium. As emphasized above 
the focus of this article is to provide a real-time description of non-equilibrium phenomena. 
Since the real time effective action within this setting has not yet received attention, we 
begin by studying the equilibrium situation first to compare with previous results in the 
literature and then continue with the focus of our program, the non-equilibrium aspects. 

In order to obtain the aforementioned effective couplings we integrate out the fermions 
to one-loop to get the (CTP) effective action for the pseudoscalars and the photons up to 
this order. Formally this is a rather straightforward procedure, the full real-time effective 
action being given by: 

expiSefflA^if^] = 

[Dtfj^V^] exp[i J d 3 x dt{i^ + {i$ - 2ga )ilj + - (+ — ► -)}] x (2.5) 
x exp[z f d 4 x{e^ + 7 / M+V> + - 2gip + (x + + ^7 5 )^ + - (H ► -)}]. 



The one-loop contribution is obtained by simply expanding the interaction part exp [iuj n t] 
to 0(e 2 g) and integrating out the fermions yielding the following expression for the ^77 
interactions: 

S n yy = —2ge 2 J d A x d 4 xi d 4 x 2 x 

x {7r + (x)A+(x 1 )At(x 2 ) Tr[S ++ (x 2 ,x) l5 S ++ (x,x 1 ) 1 < i S ++ (x 1 ,x 2 ) 1 u ] (2.6) 
-n+(x)A+(x 1 )A-{x 2 ) Tr[S- + (x 2 ,x) l5 S ++ (x,x 1 )rS + -(x u x 2 )Y} 
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-n+(x)A-^ Xl )At(x 2 ) Tr{S ++ (x 2 , x) l5 S + ~(x, x 1 )-fS- + (x li x 2 )Y] 

+n + (x)A^x 1 )A-(x 2 ) Tr[ 1 S- + (x 2 ,x) 75 5 + -(a;,x 2 )7^-(x 2 ,a; 1 )7l - (vr + -> tt")} . 
where S~ + , , S ++ and S are the real time fermionic propagators given by: 

S^ix, x') = -i (^ ± (x)^ ± (x')); (2.7) 

s- + (x, x') = -i (v>~(x)vS + (0>; s + -(x, x') = % (i) + (x')i)-(x)); (2.8) 

S ++ (—\x,x') = -S- + ( + -\x,x') Q(x -x' )-S + - { - + \x-x') Q(x' -x Q ). (2.9) 

Note that the interactions have induced cross-couplings between fields on different 
branches of the CTP contour. This means that answers to physical questions such as how 
the pion couples to classical electromagnetic sources will be affected by the interplay between 
all these cross-couplings in the effective action. In particular the equation of motion for the 



expectation value of the pion field will require all of these contributions f27]j28[ . Evaluation 
of the various terms in the effective action is more convenient in momentum space and in 
the equilibrium case we can introduce the Fourier-transformed fields and propagators: 

Ti(k) = J d 4 x e ik ' x n{x) , (2.10) 
A,{k) = J d 4 x e ik ' x A^{x) , (2.11) 
S(k) = f d*x e iHx ~ x,) S(x, x') . (2.12) 



The effective action in momentum space then takes the following form 
<i 4 pi d 4 P2 



S 



7T77 



2e*g [ f P t ° P t x (2.13) 
y J (27r) 4 (2vr) 4 V ; 



x^ + (- Pl -p 2 )A;( Pl )A+(p 2 ) J 0-Jr[S++(k- Pl ) l5 S ++ (k + p 2 )YS ++ (k)Y} + 
-^ + (-pi-P2)A+( Pl )A-(p 2 ) J 0-j T [S-^(k- Pl ) l5 S ++ (k + P2 )YS + -(k)Y} + 
-^{- Pl -p 2 )A-{p x )At{p 2 ) J 0^Tr[S ++ (k- Pl ) l5 S + -(k + P2 )YS- + (k)Y} + 
+n + (- Pl - P2 )A-( Pl )A-(p 2 ) J ^J I [S-+(k-p{) lb S + -{k + p 2 ) 1 »S--{k)Y]}. 
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^From the definitions of the real time propagators ( |2.9| ) it is clear that the four propaga- 
tors satisfy the identity: 

S ++ + S + - + S- + + S— = 0, (2.14) 

which of course continues to hold for the Fourier-transformed propagators as well. So we 
can rewrite them in terms of three independent functions in the Keldysh [^] notation, 

S R (k) = S ++ (k) + S + -(k), (2.15) 
S A (k) = S ++ (k) + S- + (k), (2.16) 

S F (k) = S ++ {k) + S—{k). (2.17) 

Here Sr, S a are the retarded and advanced propagators, while Sp is a symmetric combination 
which is closely related to the spectral density in thermal equilibrium: 

S F {k) = [1 - 2n f (\k \)]sga(k )[S R {k) - S A {k)] (2.18) 

MM) = ^ TT - (2.19) 

This relationship can be shown to arise via the Lehmann spectral representation for states 
in thermal equilibrium and it does not hold out-of-equilibrium. All these definitions and 
relations are exact and continue to hold for the full field theory in thermal equilibrium. 



B. The 'bare' ^77 vertex: 

The anomalous couplings in the CTP effective action are to be thought of in terms of 
the triangle diagram in Fig. 1 with insertions of external fields at the vertices, carrying all 
possible combinations of + and — labels on them, corresponding to the relevant portions of 
the CTP contour ^7^2^. 

Since we want to study the low-energy effective couplings, we must first clearly identify 
the scales involved in the problem. The external photons have 4-momenta pi,p2 (see Fig. 1 
for example) while the pion has momentum — p\ — j>2 by momentum conservation. Henceforth 
we will collectively denote the momenta of the external fields by -P e xt- By low-energy vertices 
we will always mean effective interactions where -P e xt << T, and since we want to remain 
close to the chiral phase transition, we choose m q « T as well. From our subsequent 
analysis it will become clear that other scales such as g T, g 2 T, originating dynamically 
from the thermal contributions to the fermion self-energy will also be extremely important. 
(Recall, in this connection our initial assumption that g << 1 so that there is indeed a 
possible hierarchy of scales T,gT.., etc.) 
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FIG. 1. The triangle diagram at finite temperature with bare fermion propagators in the in- 
ternal lines. 



First consider the regime gT << m q « T and P ex t << m q- in order to provide 
motivation for the other scales that will be important at finite T let us first calculate the 
'bare' 7177 vertex, where by 'bare' we are referring to an evaluation of the terms in the 
effective action in Eq. ( |2.6|) with undressed finite temperature propagators for the fermions 
in the loop (Fig. 2). At finite temperature the free fermion Green's functions in the real-time 
approach are: 

2ni rif(u)k)5(k 2 — m 2 ) 



S ±± {k) = (jk + m q ) 



' k' 1 



m 2 ± ie 



S T± (k) = 2m(/k + m q )5(k 2 - ml) [Q(±k Q ) - n f {uo k )} 



(2.20) 



(2.21) 



«/(M) 



1' 



k 2 + m 2 . 



(2.22) 



The local form of the low energy vertex n + A + A + (from Eq. ( |2.13p ) is obtained by 
setting -P ex t = in the computation of the triangle diagram, (as was done in the imaginary 



time formulation in |TB,n,JT|). However, in the real time formulation, setting the external 
momenta to zero inside the quark loop leads to a divergent result: 

d 4 px d 4 p 2 



5. 



7777 



im q e g 



d 4 k 



(27T) 



(2tt) 4 (2tt) 4 
1 



x {1 p laPw n + (- Pl - p 2 )A+( Pl )At(P2) 



(2.23) 



k 2 



m z + ie 



+ 2-ni nf(uj k ) 5{k 2 — m 2 



This integral has two different sources of infrared divergences. In the chiral limit m q — > 
the temperature independent term becomes sensitive to the infrared modes. At finite 
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temperature this divergence must be cured by resumming the contributions from HTL. On 
the other hand the temperature dependent terms lead to divergences ~ T 2 /e 2 . This is a 
pinch singularity that originates from the propagation of on-shell fermionic intermediate 
states present in the bath. At tree level, these quarks in the plasma propagate without 
damping over arbitrarily long times. However, in an interacting plasma the quarks will 
acquire a width T k and will only propagate during time scales ~ 1/T fc . Including this effect 
in the fermion propagator will cut off these pinch singularities. 

This naive calculation is obviously flawed but clearly indicates that using bare propagators 
and setting P ex t to zero is inconsistent with the high temperature limit. If we imagine taking 
P ex t to zero continuously, at some value the correlator will become sensitive to thermal 
corrections to the self-energy. The real part of the thermal self-energy should play the role 
of m q in the loop while the on-shell thermal width should replace the Feynman parameter 
e. 1 

Before proceeding to a more consistent description using resummed propagators, we seek 
to establish contact with the results of Jl3| and Jl7j using the real time effective action 
obtained above. We do the calculation by keeping the external momenta non-zero, and then 
take the zero momentum limit smoothly at the end of the calculation. 

Using bare fermion propagators and defining 



CO 



k 2 + 



m 



q ' 



uJ-i 



(k 



Pi) 2 + m 2 ; 



UJ2 



(k 



p 2 ) 2 + m 2 



q ■ 



(2.24) 



n 



nf(u) ; ni = n/(a;i); n 2 = n f (uj 2 ) 



(2.25) 



we find that the tt + A+A+ coupling in the effective action is given by: 



St 



7777 



-m q e g 



d A pi d i p 2 
(2tt) 4 (2tt) 4 



x e a ^ p la p 2 p n + (- Pl - P2) A+( Pl ) A+(p 2 ) x 



d 3 k 1 

(27r) 3 LOUJ1UJ2 



(1 - ni)(l - n 2 )n 



1 



+ 



— uj — p\ — ie) (0J2 — uj + p% — ie) 
1 



(2.26) 



{0J2 - uj - p\ - «e) {uJ\ - uj + pi - ie) {uj 2 - uj - P2 - «e) {uj 2 + UJ\ - p% - p\ - «e) 

1 1 

+ 



+ 



{uji - uj - p\ - ie) (uj 1 + uj 2 - Pi - P2 - ie) (ujx - uj + p± - ie) (ujx + uj 2 + p\ + p\ ~ ie) 

mn 2 (l — n) (c.c.) + 



(uj 2 - uj + P2 - ie) (out + ^2 + P1+P2- i e ) 



1 One of us (S.P.K.) would like to thank S. Jeon and L. Yaffe for pointing this out. 
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m(l - n 2 )(l - n) 



(uj - uji + pi - ie) (uj + uj 2 - P2 - ie) 



+ 



{u + U2+P2- ie) (uj — out - pi — ie) (uj 2 + uj + p\ - ie) (co> 2 - + p\ + p° _ *e) 
1 1 



+ 



(uj - <jJ\ + p\ - ie) (uj 2 - UJ\ + Pi + p\ - ie) (uj — oj\ — Pi — ie) (uj 2 - UJ\ - p\ - p® 

7 : r TT~~ m ~ (^-ni)n 2 n[c.c.] + 

(uj 2 + uj - p\ - %e)[yj2 - Wx—P\-P& - «e) 



-ie 



[1 - ni)n 2 (l - n) 



1 



(uJx + u) — p\ — ie) (uj - uj 2 + p 2 - ie) 



1 



+ 



1 



(uj — uj 2 -p 2 — ie) (uj + uji + pi - ie) (uj - uj 2 - P2 - ie) (uJ\ - uj 2 - Pi - p\ - ie) 

1 1 

+ 



(uj + UJx — p\ — ie) (uji - uj 2 - Pi - P2 - ie) (uj + uj 1 + p\ - ie) (uji - uj 2 + p\ + p\ 

- t. , . . t. t. — - ni(l — n 2 )n fee] + 

(uj -uj 2 +p 2 ) -ie)(uj 1 -uj 2 +p1+p 2 ) -ie)\ L J 

1 

nin 2 n 



-ie) 



(uj + uji + p\ + ie) (uj + uj 2 - p\ + ie) 
1 



+ 



1 



(uj + uj 2 + P2 + ie) (uj + uj 1 - pi + ie) (uj + uj 2 + p° + i e ) 6^1 +UJ2+V1+V2 + ie) 

1 1 

+ 



(uj + uji + pi + ie) (uji +U2+P1+P2 + ie) (uj + uji - p\ + ie) + u; 2 - P? - p 2 + 

- 7; — 7; 7; — (1 — ni)(l — n 2 )(l — n) Ic.c.l 1 

(uj + uj 2 -p°2 + ie)(uj 1 +uj2-p 1 ~p 2 + ie) 1 1A 2A ;L J J 



The above interaction is just one of the couplings that appears in the real time effective 
action (others are it + A~ , n + A~ A~ , etc.). An important point to note is that aside from 
the direct interactions with photons it also contains discontinuities from the on-shell propa- 
gation of fermionic intermediate states that appear in the bath and couple to the external 
electromagnetic field. These are encoded in the ie terms of the energy denominators and give 
rise to pinch singularities when external momenta are set to zero as previously discussed. 

We now make contact with the results of [|H|[17],|]1| by calculating the decay amplitude 
T(ir —>■ 77) which can be obtained from the above expression by simply ignoring the discon- 
tinuities in this effective vertex. The resulting expression can be interpreted as a particular 
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cut of the 3-loop self-energy graph |T7| (see Fig. 2) for the 7r where the photons appear as 
final states. 

In the limit of P ex t « T and P ex t « m q , the distribution functions and the energy- 
denominators can be expanded in powers of P ex t/T and P e xt/ m g respectively to yield the 
amplitude (or the imaginary part of the retarded self-energy) for 7r — > 77. After long-winded 
algebraic manipulations we obtain the following form for the vertex in the zero-momentum 
limit: 




FIG. 2: The 3-loop self-energy cut to yield the required decay amplitude in [17] 



r(vr -> 77) = m q e 2 g J 



d 3 k 1 



x 



3(1 — 2n) dn(u) k ■ pi 
2uj 2 du 2uj 

2k -p x 



(2tt) 3 w 3 

1 1 



(2.27) 



+ r- 



k-p! k-pi k-(pt+p 2 ) k-(p!+p 2 ) 
2k ■ pi 2k ■ p 2 2k ■ p 2 



k -pik ■ (pi +p 2 ) k • pxk ■ (pt + p 2 ) k ■ p 2 k ■ {p x + p 2 ) k ■ p 2 k ■ (px + p 2 ) 



dn{uj) uj 2 p{ — (k ■ px] 
du 2uj 

k-px 



k ■ px 



k-px 



+ r 



(k ■ px) 2 k ■ (px + p 2 ) (k ■ px) 2 k ■ (px + p 2 ) 
k- px k-p 2 



k ■ px(k ■ (px + p 2 )) 2 k ■ px(k ■ (px + P2)) 2 k ■ p 2 (k • (px + p 2 )) 2 
k-p 2 1 1 



k ■ p 2 (k ■ (px + p 2 )) 2 k ■ pxk ■ (px + p 2 ) k ■ pxk ■ (px + p 2 ) 



1 d 2 n{u) -> 



2 dco 2 



(k-pxY 



+ - 



k ■ pxk ■ (px + p 2 ) k-pxk- (px+p 2 )_ 



+ (pi <-> P2) 
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where uj = \jk 2 + m 2 , k — (uj, —k) . (2.28) 

The symmetrization with respect to p\ and P2 applies only to the momentum-dependent 
terms. This is one of the important results of this paper. Using this expression one can cal- 
culate the required amplitude for different external kinematical configurations. For spacelike 
photons (pi) = (P2) = and in the limit m q /T — > the expression turns out to be finite 
and is given by: 

, 9 r d 3 k 1 f 3 r 3 dn(k) d 2 n(k)) 

F(, 77 ) = ro „ e ^ / _ - - 2n( M + - -±i - -Aij . (2.29) 



This expression agrees with the result of |T7[ and after careful integrations by parts can 
be easily shown to yield the zeta-function dependent result of the Euclidean calculation. It 
behaves as ~ m q e 2 g/T 2 in agreement with [|l3l . |l~l ■ 

In the case of pion decay at rest to two back-to-back on-shell photons we find after 
straightforward substitution 

r(vr -»■ 77) = (2.30) 
d 3 k 1 ( 3 1 dn(u) d 2 n(u) (k ■ p) 2 



m « e2g te [1 - 2n( " )] 



(2tt) 3 uj 3 1 2u 21 n cu duo dtu 2 u 2 p 2 - {k ■ p)< 



This coincides with the expression found in [17] and can be shown to be equivalent to the 



results obtained in [0. In this case and also for other generic kinematical configurations 



Eq.( p.27|) in fact behaves as ~ e 2 g/T. These results are, however, valid only in the regime 
where T » m q » gT so that we can ignore the hard thermal loops. 

While this calculation seems to indicate that self-energy corrections will not play a role 
in the amplitude, the terms in the effective action also involve additional contributions 
from thermal discontinuities. These are the terms that yield pinch singularities in the zero 
momentum limit, as can be easily checked by setting the external momenta equal to zero in 
the exact expression for Eq. ( |2.2(j| ). 

The result obtained for the finite temperature vertex in the imaginary time formulation 
in pl3| , |14| , |3Tf cannot be extrapolated to the real time domain. In fll3||T4| the calculation was 



performed by setting the external Matsubara frequencies to zero. In order to obtain the real 
time behavior, the external Matsubara frequencies must be kept non-zero in the evaluation of 
the triangle diagram and the resulting expression will therefore be a function of the external 
frequencies u n defined at the discrete points along the imaginary axis in the frequency plane. 
The real-time vertex is then obtained by the continuation uj n — > i{u ± ie); the imaginary 
parts resulting from the ±ze will reveal the different discontinuities and therefore the ensuing 
pinch singularities in the limit of soft external four-momentum. 
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C. The dressed vertex with HTL resummation: 

From the above expressions it is clear that the amplitude is sensitive to momenta ~ 
m q and as m q /T is taken to zero we have to include self-energy corrections to the quark 
propagator. Furthermore we have seen that the interaction terms in the effective action also 
include the discontinuities from the quarks in the loop and these in turn will be sensitive to 
the thermal width on-shell when the external momenta are sufficiently small (i.e. P e xt << T 
where V is the on-shell width for the fermionic quasiparticles.). 




FIG. 3 . The triangle diagram at finite temperature with resummed fermion propagators in 
the loops. 

In the high temperature, long-wavelength limit, i.e. when m q ; P ex t 9^ HTL 
resummed propagators must be introduced for the quarks in the loop (see Fig. 3). The 
theory obtains HTL corrections both from the gauge bosons and the scalar sector invoving 
the a and the ir. However, since we have chosen e << g we will ignore the contributions from 
the photons to HTL. The a and the ir contribute equally to the thermal self-energy for the 
soft quarks. In addition there are also the usual zero temperature logarithmic corrections to 
the fermion mass m q which we shall ignore since these are not relevant for the discussion. The 



retarded self-energy corrections from the scalars were obtained in [j^] and we summarize the 
main ingredients that are relevant to our discussion here, in the limit m q = 0. The one loop 
retarded self-energy for the quarks is given by X(fc — «e, k) = S CT (A;o — ie, k) + Y?(ko — ie, k) 
with 



S a (A; - ie, k) = i 7o k q'(k - ie) + 7 • k e\ 1} {k Q - ie). (2.31) 
for a = a, ir. We define 



fd^-^] f. 1 ) pf(s) , 



14 



using the one-loop spectral densities given by the expressions: 



k J {2-ny i 2uj k+q uj q 

[8{s - co k+q - u q )(l + n k+q - n q ) + 5(s - u k+q + u q )(n k+q + n q )\ 

m / x , 2 f <Pq 1 k ■ q 



k J {2^2uj k+q uj q k 2 

[S(s - io k+q - u) q )(l + n k+q -n q )-S(s- u k+q + Cu q )(n k+q + n q )} , (2.33) 



where n^ + ^ , n q are the Bose-Einstein and Fermi-Dirac distribution functions and 00 , 00 are 
the bosonic and fermionic frequencies respectively. The Bose-Einstein distributions for the a 
and 7r fields are functions of their respective masses. We have also neglected the ultraviolet 
counterterms as these are irrelevant for our discussion. 

The terms proportional to S(s — u k+q — uj q ) in the spectral densities correspond to the 
processes q <-> + q, where is either aaorir field, and determine the usual two particle 
cut that survives in the zero temperature limit. The terms proportional to S(s — u) k+q + u) q ) 
correspond to Landau damping and processes of the type <-> g + g, which can only occur 
in the medium. Near the critical temperature and in the strict chiral limit, m a ,m n — > as 
T — > T~ and m q — > 0. Therefore, in this limit, the Landau damping terms lead to a HTL 
contribution to the quark propagator. 

We extract the HTL contribution to the quark self-energy in the strict chiral limit and 
near the critical temperature so that T >> m q ,m cr ,m 7T , P ex ^ in the usual manner |32| |36 
and find 



Z HTL (k - ie; k) 



M eff 



7 ln 



k + k 



16 



+ 7 • k 



2 _ h. i n k ° + k ~ ie 
k ko — k — ie 
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2 rp2 



(2.34) 



This HTL contribution is very similar to that of the gauge fields. In fact including the gauge 
field is tantamount to replacing g 2 — > g 2 + e 2 in the above expression. 

In the HTL approximation, the effective quark propagator has poles describing two 
branches of collective modes In this approximation the self energy only has an 



absorptive part below the light cone, so that to leading order in the HTL limit the on-shell 
width of the collective modes vanishes. The effective mass for the collective modes M 2 ^ does 



not break chiral symmetry |3"2"|-|3T)|] but serves as an infrared regulator. In order to obtain 
the width of the quasiparticles to cure the pinch singularities we must go beyond the HTL 
contributions. 

Although including the damping of excitations in the medium and consequently a width 
in their spectral densities is a physically reasonable manner to cure the pinch singularities, 
an alternative approach has been discussed in the literature |37 . 
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D. Including the width - beyond HTL: 



There are several different contributions to the thermal width of the quark in the medium. 
These come from scalars, pseudoscalars and gauge bosons and have very different physical 
origins. We study separately the contributions from the scalar and pseudoscalars and of 
gauge bosons to clarify the differences. 



1. Scalar /pseudoscalar contributions: 



In the expression for the one loop self-energy given above ( |2.31| )- p.33| ) the terms propor 



tional to S(s — ujk+q + &q) could lead to an imaginary part on the mass shell of the collective 
modes. This would be the case if the delta function is satisfied for s = u>±(k) with u±(k) 
being the dispersion relation for the different branches of collective modes. This relation 
clearly depends on the mass of the scalar or pseudoscalar particle in the loop. Since this 
delta function arises from the process <-> q + q the constraint is satisfied if the scalar (or 
pseudoscalar) boson can decay into a qq pair. 

The fact that the decay of the scalar leads to a thermal width for the fermion has been 
previously recognized in p2 and can be understood from a simple kinetic argument: the 



change in time of the spin-averaged quark distribution can be written as a balance equation 
with a gain term minus a loss term. The gain term arises at lowest order from the decay of 
the scalar into qq pairs with a factor 

\M fi \ 2 n k+q (l - n k )(l - n q ), 

while the loss term arises from the 'recombination' of quark-antiquark pairs into a sigma 
meson with probability 

\M fi \ 2 {l + n k+q )n k n q 

with |A^/j| 2 the usual transition matrix element. The linearization of the resulting Boltz- 



mann equation leads to twice the damping rate p2| of the quark of momentum k. Thus 



the decay of the scalar into qq pairs in the medium induces a width for the quark or the 



collective excitations 22 



In the strict chiral limit with T «T C , the mass of the a meson « 660 MeV while that of 
the constituent quark is ~ 300 MeV so that the kinematical decay condition could be fulfilled. 
Near the critical temperature the sigma mass vanishes as does the constituent quark mass 
(in the strict chiral limit) but now for T » m q the important mass scale is determined by 
M eff which for g « 1.5 , T « T c « 150 MeV yields M eff « 80 MeV. Whether the constraint 
could be fulfilled requires a detailed study of the masses and their temperature corrections 
for the quarks and the scalar mesons, a task beyond the scope of this article. However for 
T ~ T c where in the strict chiral limit m q = and M e // ~ 80 — 100 MeV, the sigma meson 
will be able to decay into collective excitations giving the quarks a thermal width of origin 
very different from the more familiar collisional width. From the one-loop expression for the 
scalar and pseudoscalar contributions to the quark self-energy we see that such a width will 
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be of order g 2 T (for a detailed expression for the width see [0])- In the perturbative analysis 
this width is much smaller than the effective mass of the collective excitations ~ gT. 

If the temperature dependent mass of the sigma prevents its decay into collective modes 
by kinematics, a collisional contribution to the width of the quark appears at two loops in the 
quark self-energy. The argument given above is valid whenever the fermion momentum in 
the self-energy loop is » gT. In particular, for T < T c but for m a (T) >> gT, the on-shell 
delta function that determines the width requires that the momentum of the fermion in the 
self-energy loop be > gT [22] and the use of the bare quark propagator is justified. On the 
other hand very near the critical temperature whenever m cr (T) « gT for soft external quark 
momentum the internal quark line will also be soft in the kinematic region that contributes 
to the width. In this case a full HTL resummed quark internal line must be used. The 
analysis of such case has been provided in detail by Thoma [21] who concluded that the 
width of the quark in this limit is also given by ~ g 2 T. This result includes the collisional 
width but in the HTL approximation for the intermediate fermion in the self-energy diagram 







In summary, for scalars (and pseudoscalars) the thermal width of the quark collective 
modes has two different contributions both of order g 2 T: the decay of the scalar (sigma) 
into pairs of collective modes if kinematically allowed and collisional damping in the HTL 
limit . 



2. Gauge boson contribution: 



The fermion damping rate arising from gauge boson exchange has been studied in detail in 
the HTL approximation by several authors P8"|-f4"0|] . The gauge contribution to the quark self- 
energy presents infrared divergences in bare perturbation theory. The infrared divergences 
associated with the exchange of longitudinal gauge bosons are a result of small angle Coulomb 
scattering and at zero temperature are those of Rutherford scattering. At finite temperature 
the longitudinal gauge boson (instantaneous Coulomb interaction) is screened by the Debye 
screening mass mo oc e T which cuts off the infrared and leads to a finite contribution to the 
damping rate from longitudinal gauge bosons. For soft external momentum this contribution 
has been found to be given by [ p^ , p^ j 



r l ' 9 (k) — a AT , (2.35) 

= e 2 /4vr. 

For a fermion excitation at rest the damping rate has been computed in [38,41,42] and 
found to be given by 



with A a constant that can be found numerically |38| and in QCD is of 0(1) and 



T 9 (k = 0) = a C T 
with C again being a numerical constant of 0(1 



(2.36) 



The detailed analysis of Refs. [ 3S , 41 



@E1- 

in QCD lead to the result that the damping rate 
of moving quasiparticles with momenta k 3> e 2 T in a non-abelian plasma are given by (up 
to an overall constant that depends on the gauge group structure) 
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T 9 ± {k) = ^\v ± {k)\ log- , (2.37) 
Air e 

where v±(k) are the group velocities for the two different branches of collective modes. This 
expression is not valid for k = where the damping rates of the fermion and plasmino at 
rest do not have the logarithmic behavior in terms of the gauge coupling |38],f|l], nor near 
k = k m i n where the group velocity of the plasmino branch vanishes. Since the HTL structure 
in QED is similar (up to gauge group factors) to that of QCD (and scalar QED), Pisarski 
p8| suggested that a similar form of the damping rate should be valid in a QED plasma, 
despite the fact that there is no magnetic screening in the Abelian theory. 

In QED the transverse photon propagator is only dynamically screened via Landau damp- 
ing and the infrared divergences remain, possibly to all orders in perturbation theory. More 
recently [^] a detailed study of the fermion propagator in the Bloch-Nordsieck (eikonal) 



approximation in real time revealed that for lud t \v±(k) \ 1 

S k (t) « e - aTlv±{k)l tln(uJD 1 lv±m , (2.38) 

with ujd ~ eT the Debye frequency and v±(k) the group velocity of the fermionic collective 
modes. Although this is not an exponential relaxation that would emerge from a Breit- 
Wigner resonance shape of the spectral density, it does reveal a particular time scale from 
which a damping rate can be extracted and is given by 

r±(/c) «aT |f±(fc)| In - . (2.39) 
This result has been recently found in scalar QED (which has the same HTL structure as 



QED and QCD to lowest order) within the dynamical renormalization group method |4(| 



At this stage it is important to describe the physics that leads to the damping rate from 



the gauge boson contribution pSj-pOl, which is rather different from that of scalars found 
before. The constraints from energy momentum conservation in the imaginary part of the 
self-energy can only be satisfied below the light cone where the HTL resummed gauge boson 
propagator has support that arises from the Landau damping cut. The infrared process that 
leads to the non-exponential relaxation is the emission and absorption of soft photons at 
almost right angles with the moving fermion p8| , p9| . 

Thus we summarize the contribution to the damping rate of the collective excitations by 
the HTL resummed gauge boson exchange: 

r 9 Jk) « a T \v ± (k)\ hi- for k^>g 2 T, (2.40) 

e 

r 9 ±(k) ^aT for k = . (2.41) 

The HTL resummed retarded and advanced propagators for the soft quarks can be sum- 
marized as follows (here we also include the gauge contribution to the HTL self-energy) 
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S R {k ,k) 
S A (k ,k) 



7° A - 7 • Jfe 5 
7° A* - 7 • fc 5* 



.4 



*2 



*2 



(^ 2 + e 2 )T 2 
8k 



. k ±k k + k-ie 
=Fl ± — — — In 



2A- 



k n — k — ie 



+ U ± (k ,k) . 



(2.42) 

(2.43) 
(2.44) 



where n±(/co, k) are the (9 (g 2 T, aT, aTln corrections to the HTL self-energy that deter- 
mine the width of the quark. The retarded and advanced propagators uniquely determine 
the time-ordered propagator via the spectral representation: 



S++ = S R + S A + [1 - 2*1/(1**,!)] sgn(A;o) (S R - S A ) 



(2.45) 



The anomalous vertex with dressed fermion propagators is shown in Fig. 3. Using the 
general expression fl2.13|) , incorporating the on-shell width ~ g 2 T into our definition for A 
and assuming that the external momenta are negligible compared to all other scales in the 
problem we find that the n A+ A+ vertex is given by: 



C+ + + ~ 

7T77 ~ 



d A k 



(27T) 



64m„ e g 



Re- 



2 j d^pi d A p 2 



(2tt) 4 (2tt) 4 
1 



x {t Pla p 2p n + (- Pl - pa) i;( Pl ) A+(p 2 ) (2.46) 

1 



A*{k ,k)-B%k ,k) 



+ [1 -2n / (|fco|)]sgn(fc )Im 



A 2 (k ,k) — B 2 (k ,k) 



The infrared region for the real part is cut-off at the scale (g 2 + e 2 )T where we have 
included the HTL contribution from the gauge fields, and the imaginary part now has support 
below the light cone arising from the HTL and also near the mass shell of the collective modes 
because of the width T w (g 2 T, aT, aTln -) which now cuts off the pinch singularity. The 
real time vertex is now infrared finite and free of pinch singularities and the local limit of 
zero frequency-momentum in the quark loop can be taken safely. 

At this point we must make an important comment - in considering contributions be- 
yond the leading order in the HTL program, we must also include consistently the vertex 
corrections. There are in fact two sources of vertex corrections coming from a, n and photon 
exchange. As long as g » e, it is the former that is more important. The inclusion of the 
photon exchange is necessary as a matter of principle to ensure that the Ward identities are 
satisfied and gauge invariance is respected. Including the vertex correction along with the 
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quasiparticle width makes a detailed evaluation of the effective vertex an extremely difficult 
task, certainly beyond the scope of this article. The conclusions that we will obtain are 
therefore tentative pending inclusion of the vertex corrections. 

As mentioned in the introduction one of our principal goals is to understand the effective 
interaction of the neutral pion with gauge fields to assess the possibility of enhancement 
of neutral pion condensates in the presence of classical electromagnetic fields produced in 
peripheral collisions. To obtain the interaction of the pseudoscalar with classical background 
fields we simply replace the fluctuating fields in Eq. ( |2.13|) with c-number expectation values 
or classical sources. This yields the following coupling of the pion to classical fields A^, 



= J |^|^ + (-^ - A M MV*) J ^ x (2.47) 

{Ti\s ++ (k - Pl h 5 s ++ (k + P2 h ll s ++ (k)Y} 

-Tr[S-+(k - Pl h 5 S ++ (k + p 2 h> 1 S + -(k)Y} 
-Ti[S ++ (k - Pl )^S+~(k + P2 )YS- + (k)Y] 

+Tr[ 1 S- + (£; - Pl ) l5 S + '(k + p 2 )Y S~~ (k)Y}} 

Thus the interaction of the pion with classical gauge fields is determined by a combination 
of all the different terms in the CTP effective action. Rewriting this expression in terms of 
the propagators in the Keldysh notation |eqs.( |2.15| )-( |2.17| )1, we obtain 



S nll = -2e 2 g I ?r + (-pi - p 2 )A tl (pi)A v (p 2 ) 



d A k 



(27r) 4 (27r) 4 v ^ uy ^ J J (2tt) 4 

{Tr[S R (k - Pl )j 5 S R (k + p 2 )^S R {k)Y] + Tr[S A (k - Pl ) l5 S A (k + p 2 )^S A (k)j"} + 



(2.48) 



+Tr[S A (k - Pl )j 5 S F (k + p 2 )^S A {k)Y) + Tr [S A {k - Pl )>y 5 S R (k + p 2 )rS F (k)Y} + 



+Ti[s F (k - Pl ) l5 s R (k + p 2 )rs R {k)YW 



The local limit of the pion-photon vertex is obtained by taking the external frequency 
and momentum to vanish in the quark loop. After the HTL resummation and accounting 
for the width of the intermediate quarks, this limit is unambiguous and finite. Thus we find 

<? — 

16m, e 2 , / x ji^ p la p w *+(- Pl - P2 )i+( Pl )i+(p 2 ) / y^I* 
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' 1 [l-2n/(|*b|)]sgn(fco 



(A 2 ~B 2 ) 3 (A* 2 -B 



*2\3 



(A 2 - B 2 f (A* 2 - B* 2 f 



+vertex corrections (2.49) 

From this result the strength of the local pion-photon vertex can be extracted, which is 
free of ambiguities and non-singular. It is of the form 



2 T 

5 , 7T77 rp2 6 9 J~ 



g,a,alog-, 
e 



(2.50) 



The function T is finite and dimensionless and certainly quite complicated by the loga- 
rithmic dependence on the coupling of the width and by the vertex correction. It is beyond 
the scope of this article to obtain a closed or approximate expression for this function. The 
main point of writing the eqn. ( [2.50|) is to emphasize that after curing the pinch singularities 
by including the width (and vertex corrections) the vertex vanishes at T = T c in the strict 
chiral limit. This is in agreement with the results in [DJ. The necessity of including HTL 



corrections was pointed out by In the present context, the CTP effective action requires 
the inclusion of a width as well, in order to regulate the pinch singularities. This is one of 
the important results of the present work. 

A crude, order of magnitude estimate for the vertex may be found in the narrow-width 
approximation, assuming that the width T ~ 0(g 2 T). Then, the dominant temperature 
dependent term in g nll is ~ e 2 gm q x (T 2 /r 4 ) which is 0(e 2 /(g 7 T 2 )). (The T 2 /(T 4 ) factor 
may be easily understood by going back to the discussion after Eq. (2.23).) 

Before we embark on the study of nonequilibrium effects it is convenient to summarize 
the conclusions for the equilibrium case 

• The real time interaction vertices for the pseudoscalar and photons obtained by inte- 
grating out the quarks in the real time formulation display infrared divergences in the 
local limit when the bare quark propagators are used with m q = 0. These divergences 
are cured by an HTL resummation of the quark lines and by considering contributions 
beyond HTL that give an on-shell width to the collective excitations. The resummed 
quark propagators now give rise to an unambiguous and finite local interaction vertex. 

• The computation of the effective vertex in imaginary time in terms of bare propagators 
[|ll]|,|n|,[n[] cannot be used to discuss either the decay amplitude of the neutral pion or 



any dynamical aspects. In that computation the external Matsubara frequencies were 
set to zero and therefore the analytic continuation necessary to obtain the real time 
vertices cannot be performed. This computation cannot reveal the subtle infrared and 
pinch singularities that arise from finite temperature contributions to the discontinu- 
ities of the triangle diagram that are responsible for ambiguities in the limit of soft 
external momenta. A resummation of the quark propagator including both HTL con- 
tributions and the corrections that lead to the width of the fermionic collective modes 
is necessary to regulate the infrared divergences in the soft limit. 
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The infrared and pinch singularities of the real time vertex imply that the local limit 
of the pseudoscalar-photon vertex is very sensitive to the details of long- wavelength 
physics and is therefore non-universal in the sense that it depends on the details of the 
underlying model. This is in striking contrast with the situation in the vacuum where 
the vertex in the local limit is completely determined by the triangle anomaly and 



the Adler-Bardeen 'theorem' [43| guarantees that the vertex does not receive radiative 



corrections. At finite temperature the situation is more complicated. As originally ob- 



served in [II |, while the anomaly equation is independent of temperature and does not 



receive radiative corrections, the vertex itself is not determined solely by the anomaly. 



Our studies in equilibrium complement those of Gelis |T7[ and lead to the conclusion 
that the vertex is sensitive to the infrared and thus to the details of the low energy 
scalar and pseudoscalar sectors. 

After the HTL resummation and accounting for higher order corrections that determine 
the fermion width the local limit of the vertex is given by Eq. ( |2.50| ) which vanishes in 



the strict chiral limit m q = in the symmetric phase, confirming the results in |]l3| , p 
Hence, despite the fact that the imaginary time computation of the vertex had missed 
the infrared sensitivity, it does reveal the correct behavior for the pseudoscalar-photon 
vertex in the strict chiral limit near the phase transition, when the constituent quark 
mass vanishes. 

The full benefits of the real-time effective action will be appreciated in a truly out of equi- 
librium situation. In the earlier sections we reproduced previous equilibrium results within 
this approach and pointed out certain important subtleties pertaining to the appearance of 
pinch singularities. We now turn our attention to the true non-equilibrium situation which 
is the focus of our work. 



III. 7T77 OUT OF EQUILIBRIUM 

We expect that nonequilibrium phenomena will be the rule rather than the exception in 
heavy ion collisions. In particular, in terms of DCC formation, we might expect that since 
electromagnetic effects break isospin, they might be able to bias the neutral to charged pion 
ratios in a way that could be detected in an event by event analysis. 

The scenarios of formation of DCC-type configurations envisage a strongly 'quenched' 
phase transition, in which the cooling of the plasma occurs via hydrodynamic expansion on 
a time scale far shorter than the equilibration scale of long- wavelength fluctuations [0-0. 
Within the setting of the linear sigma model this situation is modeled by suddenly changing 
the form of the potential for the scalars (and pseudoscalars) from one with an unbroken 
symmetry minimum to another that displays symmetry breaking minima. 

While most of the previous studies focussed on an effective theory of scalars and pseu- 
doscalars (mainly the linear sigma model), in order to study the effects of the anomalous 
vertex on possible enhancements of neutral pion condensates (constituent) quarks must be 
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added to the model. Following phenomenologically motivated descriptions [0,0 we adopt 
the U{2) ® U (2) gauged constituent quark model. This is, of course, a simple extension of 
the model of the last section, with an SU(2)l x SU{2)r global symmetry with two quark 
flavours and N c = 3 colors, coupled to electromagnetism: 

£ = z$(£ - ieQ4)iP - 2gi>(ar° + in ■ r 75 )</> - -F^F» V + C n . (3.1) 

Here (t°,t) are the isospin generators proportional to the Pauli matrices, satisfying the 
normalization conditions Tr(r a r 6 ) = 5 ab /2, r° = 1/2 and 

-C) 

where we have suppressed all color indices. The U(l) em charge is given by Q = r°/3 + r 3 . 
£ 7r determines the dynamics of the pions and the sigma field and is a gauged 0(4) scalar 
theory given by, 

At = \{d,af + \{d^f + ^(vr 2 + 7T 2 2 )A^ + j^ - A(a 2 + tt 2 - v Q f + ha (3.2) 



3n = e (7Ti<9 M 7r 2 - 7T 2 9 M 7ri) (3.3) 

Our approach to the problem will be to consider the evolution in time of the chiral 
order parameter/condensate (a), (where the expectation value is taken in the time evolving 
density matrix), and to study how this evolution feeds back into the fermion propagators in 
the triangle diagram. 

We consider the initial state to be described by free fields in thermal equilibrium at a 
temperature Tj > T c and evolve the equations of motion with the T = masses. This 



'quench' scenario was originally described in |qj27|. Several authors studied the cooling by 
hydrodynamic expansion |||7||| with similar results. Therefore we consider preparing an 
initial density matrix that describes the system in equilibrium at an initial temperature 
above the critical at an initial time U which for convenience will eventually be taken to be 
ti = 0. This initial density matrix is now evolved in time with a potential corresponding to 
the zero temperature situation [^7J . 

At time t = ti the temperature suddenly falls below the critical value and the potential 
develops two (non-degenerate) minima. The condensate now evolves in time from the initial 
value towards the final equilibrium point with the minimum (free) energy. We always begin 
with an equilibrium initial state where the chiral order parameter has a non-zero expectation 
value (cr(x)) oc h which is the explicit symmetry breaking term and yields a non-zero bare 
mass for the quarks. 

We should emphasize at the outset that the infrared and pinch singularities that plagued 
the equilibrium calculation are absent in this scenario. The reasons for this are 
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We begin with a non-zero condensate value which yields a finite constituent quark mass 
at Ti > T c . This arises from the explicit symmetry breaking in the linear sigma model, 
which is responsible for the physical masses of the pions. 

The range of integration time is finite, from an initial time ti = up to the time t of 
interest when the phase transition is complete, t « 10 fm/c. The fact that the time 
interval is finite introduces a natural infrared cutoff and prevents the pinch singularities 
associated with the propagation of on-shell intermediate states over very long times 



as emphasized in the non - equilibrium situation [44 . That real time acts as an 



infrared cutoff for the propagation of intermediate states has been understood within 



the framework of the damping of collective excitations in a plasma in [RS , 40 ] 



Since the scalar couples to quarks and pions, the rolling of {cr(x)) = o"o(t) has two impor- 
tant effects: i) an effective time dependent mass for the pion fields that leads to instabilities 
in the spinodal region of the potential and consequent growth of pion fluctuations ||6|-P,p7|, 
ii) a time dependent constituent mass for the quarks. 

These nonequilibrium effects result in the production of a large number of pions via 
spinodal growth of pion fluctuations p|-BL|27|,E5l and production of qq pairs through the time 



dependence of their mass [pq] . The quantum backreaction of quarks and pions affects the 
dynamics of (cr(x)) so that we cannot simply focus on the classical equations of motion and 
a self-consistent treatment of the dynamics of the condensate incorporating quantum effects 
is essential. Detailed investigations of such dynamical situations in the O(N) model and 
the scalar sector of the linear sigma model using nonperturbative schemes such as the large- 
N limit have been carried out in previous works by various groups ||7|,|9|, [27| , |29| without the 



inclusion of quarks. The out of equilibrium evolution including fermion fields is treated in 



As before, the full quantum field a ± (x) can be split up into its classical, c-number ex- 
pectation value and quantum fluctuations about such expectation value 

a ± (x,t)=a (t)+ X ± (x,t) (3.4) 
(a ± (x,t)) = a (t) , ( X ± ) = 0. 

The condensate will be introduced as a time-dependent mass term for the quarks and its 
dynamics will be treated non-perturbatively. We split the fermionic part of the Lagrangian 
into a free part and interaction terms as follows: 

£f — + £>int 

C = - ga (t)}^ + -(+ — -) (3.5) 

dm = e^ + YQA+*P + - igr ■ vr+^ + 75^ + -(+—►-) + •••■ 

where the dots stand for terms which are not relevant to the pseudoscalar-photon vertex to 
one loop in the triangle diagram. We now combine the results of the (mean-field) large N 
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equations for the linear sigma model [||[?||| |27 , 2j| with those obtained for the Yukawa theory 
out of equilibrium |28] to obtain the self-consistent equations of motion in the mean field 
approximation. 

a (t) - 4A v 2 cr (t) + 4A o*(t) + 4A cr (t) (tt 2 (x)> - g(^) = h 



^)> = 1/(08 l^(0! 2 COth 



^fc,7r(0) 

2T, 



**(*) + uf^it) Mt) = ; uljt) = k 2 - 4A ^ 2 + 4A cr 2 (t) + 4A (7? 2 (*)> 



[^-J(7o(tMl)=0 



(3.6) 



where we have set the number of pions N = 3. The reader is referred to [|^,[7tp|,p7]-p9|1 for 
details on these equations and their initial conditions. 

In the above equations we have neglected the contribution from the classical gauge fields 
to the back-reaction problem for the dynamics of the condensate. This is justified since 
the electromagnetic effects are of C(e 2 ) but A, g » e and A ~ g |f|,f7|,|9"|, |2"7| , p9| . This set of 
coupled self-consistent equations determine the time evolution of the condensate including 
the large contributions from the pion and fermionic fluctuations. Although the integrals 
associated with (n 2 ) ; (ipip) are ultraviolet divergent and renormalization would be required, 
the model under study is only useful as a low energy effective description with a natural 
cutoff A < 1 GeV. 

The fermion propagators in the background of the condensate for a 'quench' scenario 
from an initial state with temperature Tj are given by HB 



St-(t,t') = -Sf(t,t') = -i / d 3 x e"^ (^-(6,t')^ + (x,t)) 



Sz+(t,t') = SUt,t') = i d 3 x e~^ (VT(f,f)f + (ff.O> 



(3.7) 
(3.8) 



st + (t,?) = sut,t')e(t - f ) + sf(t,t')e(t' - 1) 



Sr&t) = SUt,t')Q(t' -t) + SS(t,t')0(t - t>) 



(3.9) 
(3.10) 



Sf(t,t') = S$t,t) [1 - n f M) - S<^t,t')n f (u k ) 
S$(t,t') = S<~(t,t') [1 - rifM] - S>fi,t')n f (u; k ) 



i 



(3.11) 
(3.12) 

(3.13) 
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We have listed the nonequilibrium fermion propagators for the case where the initial state 
is thermal with temperature Tj = 1/(3 and the chiral condensate has the initial value cr(0) 
for all times t < U. and S^j. are the Wightman functions for the zero temperature case 
and are given by 

Sok(ti t ) = 

-ifk(t)ti{t') [Wk(t)j -y-k + ga (t) 
Sok(ti t ) = 

-if*(t)f k (t') \wm l0 - 7 • £ - ga (t) 



1 +7o 



1 - To 



^(070-7-^ + ^0(0 
>V fc (07o-7-^-^o(0 



m(t) 
d 2 



hit) 



dt 2 
fk(t < t 



' h{t) ' 

+ k 2 + g 2 al(t) -ig& (t) 



fk(t) = 



0) 



2u k (u k + ga (ti = 0)) 



(3.14) 
(3.15) 

(3.16) 

(3.17) 



where the mode functions fk{t) determine the solution of the time dependent Dirac equation. 
Details of the solution to the Dirac equation and the derivation of the nonequilibrium prop- 
agators for the fermions are provided in Appendix B. It is this dressing of the fermions by 
the time-dependent condensate which will ultimately lead to nonequilibrium modifications 
of the 7r°77 couplings. In terms of these propagators we find [^8[ 



d 3 k 
(2+p 



Tr S<(t;t) 



2NfN c 



d 3 k 
(2+p 



l-2k 2 \f k (t)f 



(3.18) 



with Nf, N c the number of flavors (2) and of colors (3) respectively. As mentioned in the 
introduction, it is not our goal here to fully study the coupled set of equations, but to focus on 
how this mean field evolution provides nonequilibrium effects that modify the pseudoscalar- 
photon coupling. 



A. Integrating out the fermions: 

As in the equilibrium case, the effective couplings will be induced once the fermions 
are integrated out at the one-loop level via the triangle diagram. As before, we extract 
the nonequilibrium couplings by expanding exp[iSi n t] and integrating out the fermions at 
one-loop. This procedure automatically incorporates all the effects of the medium or the 
condensate which couples to the fermions in the loops and we find the following expression 
for the 7r°77 interactions in the nonequilibrium effective action: 
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S%?* = - 9 -^ J d 4 x d A x x d 4 x 2 x 

x{n 0+ (x)A+ l (x 1 )A+(x 2 )Tr[S ++ (x 2 ,x) l5 S ++ (x,x 1 )rS ++ (x 1 ,x 2 )Y] (3.19) 
-n 0+ (x)A;(x 1 )A-(x 2 )Ti[S > (x 2 ,x) l5 S ++ (x,x 1 )YS < (x 1 ,x 2 )Y] 
-7r 0+ (x)A^(x 1 )A+{x 2 )Tr[S ++ (x 2 , x) l5 S<(x, x 2 ) 1 v S > (x 2 , x x )^] 
+n 0+ (x)A^(x 1 )A-(x 2 )Tr[S > (x 2 ,x) l5 S < (x,x 2 )rS--(x 2 ,x 1 )r}}- 

where S > , S < , S ++ are the nonequilibrium fermionic propagators given by ( |3.7| - |3.17D (see 
Appendix B for details). 




G(t) 



FIG. 4: The ^77 vertex out of equilibrium. 

Fig. 4, illustrates the nature of the anomalous couplings which are induced out of equi- 
librium. Note that cross-couplings between fields living on different branches of the CTP 
contour have been induced in Eq. (|3.19| ). We note that if the zero temperature, equilibrium 
form for the fermion propagators is used the cross-couplings vanish and only the n 0+ A + A + 
vertex remains. Now the effective interactions in Eq. (|3.19|) are manifestly non-local. In 
practice a local vertex is obtained via a derivative expansion or an expansion in powers of 
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the external 4-momenta where each derivative is suppressed by the mass of the fermions 
in the loop which are being decoupled. In the usual equilibrium case the fermion mass is 
simply gvo ~ gf n and the derivative expansion is an expansion in powers of (d^/gf^) or an 
expansion in powers of (d^/T) in the high temperature regime. 

In the strongly out-of-equilibrium scenario under study, such a derivative expansion is 
unwarranted. Furthermore, the lack of time translational invariance due to the rolling of 
the condensate out of equilibrium prevents Fourier transforms in time and a subsequent 
expansion in frequencies. The Green's functions and consequently the kernels are no longer 
functions of time differences, but in fact functions of the individual time arguments. It is 
also clear that separating the center of mass and relative time variables as in the case of a 
Wigner representation is of little use since there is no obvious separation of time scales. 

However if the electromagnetic field can be treated classically and if it is slowly varying 
on the time scale of evolution of the condensate (5 — 10 fm/c ||[7||||27|]) we can extract a 
local interaction term of the type ~ fl^rrC^) tt° E ■ B where the time-dependent coefficient 
determines the strength of the anomalous coupling during the stage of nonequilibrium evolu- 
tion of the condensate. The restriction to static classical electromagnetic fields is consistent 
with the assumptions used by |20| in their estimate of anomalous enhancement of neutral 
pion condensates. 

B. Coupling to quasi-static E and B fields: 

Although the non-localities in the nonequilibrium vertex are difficult to handle in general, 
couplings to quasi-static electromagnetic fields are relatively easy to obtain and we will now 
outline the derivation of these couplings. The assumption of quasi-static electromagnetic 
fields entails that their time variation is on scales much longer than those of the nonequilib- 
rium kernels (~ 5 — lOfm/c |6||/],[|,[27[] ) that enter in the vertex. Under this assumption we 
can take these electromagnetic fields outside of the time integral. We first introduce spatially 
Fourier-transformed fields: 

Tr 0+ (q,t) = J d 3 x e~ iq ^ 7r 0+ {x,t) , 

t)= J d 3 x e-* 3 A±(x, t) . (3.20) 

In terms of these fields, we find that the following anomalous interactions are induced at 
one-loop in the nonequilibrium action: 



7f 0+ (gi,t)Aj(g 2 ,ti)A+(-gi -q 2 ;t 2 ) J —-r 3 Vl l 't + (q 1 + q2 + k,k + q 2 ,t,t } t 1 ,t 2 ) + 
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— — g^_(gi + <f 2 + fc, k + q 2 , k; t, h, t 2 ) + 



(2nf 
d 3 k 



/UK 
—— z V^_ + {q l + q 2 + k,k + q 2 , k; t, t u t 2 ) + 



j—^V^.{qi + q 2 + k,k + q 2 , k; t, t u t 2 

where, V+++, V++_, V+1 + , and are given by the following expressions in terms of the 

nonequilibrium propagators, 



Vr ++ (h,k 2 ,k) = (3.22) 
Tr[5>(t 2 ,t)75^ > 2 (^ii)7^ fe < (ti,t2)7ie(t 2 -t)0(t-ti)0(t 2 -ti) + 
Tr[5<(t 2 ,t) 7 5^ < 2 (t,ti)7^ fe > (ii^ 2 )7l0(t - fa)e(ti - t)e(ti - f 2 ) + 
Tr[5>(t 2 ,t)75^ < 2 (t,ti)7^ fe > (^i^2)7l0(t2 -t)6(ti - t)0(ti - t 2 ) + 
Tr[5>(t 2 ,t)75^ < 2 (t,ti)7^ fc < (^i^2)7l0(t2 -0©(*i -*i) + 

Tr[5<(t 2 ,t)75^ > 2 (^ii)7^ fc > (^i^2)7l0(t " * 2 )e(f - tOe^x - f 2 ) + 
Tr[5<(t 2 ,t) 75 ^ > 2 (t,t 1 ) 7 ^ fe < (ti,t 2 )7l0(t - i 2 )0(i - h)e(t 2 - h) , 

Vr + -(ki,k 2 ,k)= (3.23) 
Tr[5>(t 2 ,t)75^ > 2 (i^i)7^ fe < (ti,t2)7l0(t-ti) + 
Tr[5>(t 2 ,t)75^ < 2 (t,ti)7^ fe < (^i^2)7l0(ti -*) , 

Vr- + (ki,k 2 ,k)= (3.24) 

Tr[5>(t 2 ,t)75^ < 2 (t,ti)7^ fe > (^i^2)7l0(t2 -t) + 
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Tr[5< i (t 2 ,t) 75 S<(t,t 1 )7^>(t 1 ,t 2 )7ie(t-t 2 ) ) 

Vr—(ki,h,k)= (3.25) 
Tr[5>(t 2 ,t) 7 55<(t,t 1 )7^>(t 1 ,t 2 )7ie(t 2 + 

Tr[5>(t 2 ,t) 75 5<(t,t 1 ) 7 ^<(ti,t 2 )7l0(ti -t 2 ) , 

and ki = k + q x + <f 2 , £ 2 = A; + g 2 . 

We now assume that the photons can be treated as classical quasi-static background fields 
in accord with our original goal of studying the response of the neutral pion field to large 
electromagnetic fields in peripheral collisions and in agreement with the treatment of 



We then perform a spatial derivative expansion to extract a local coupling. The result will 
be gauge-invariant under restricted (static) gauge transformations. While such couplings 
cannot be used to obtain the pion width these are relevant for addressing the possibility 
of neutral pion condensates induced by the electromagnetic fields produced in peripheral 
collisions or through the phase transition. 

We introduce c-number expectation values (mean fields) for the electromagnetic fields, 

A das S ical = ^ = A ^ ( 3 _ 26 ) 

In the presence of these classical sources the coupling of 7r° to the electromagnetic fields 
takes the following form, 



jf* I (2^3 (2^3 ^° + (^) M&) >< (3-27) 

oo poo f rj^h 

' -I*V TrfJ,V TrflV \ 



This is obtained from Eq. ( 3.21j) by replacing the fluctuating gauge fields with their 



c-number expectation values. The range of the time integrals in the expressions above is 
from ti = to t as determined by the step-functions (see below). 

In this nonequilibrium scenario, the anomalous coupling to classical fields is then 

e 2 g^ Jdtj f? 3 * 0+ (qi,t) A,(q 2 ) A,(-qt - q 2 ) X V"fa,fct). (3.28) 
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where 

poo roc f fj^h 

V^(qnq 2 ;t)= jf_ dt x jf_ rft 2 y — 
[" {Tr [S> (t 2 , t) 75l <?> (*, ^1)7^(^1, * 2 )7 



Tr [S< (t 2 , t) 75 5< (i, tl )YS> (h, t 2 )Y] } 



xG(t - *i)e(t - * 2 )+ 

e(t - ti)e(t - t 2 )e(f 2 - *0+ 



X 



+ {n[S£(M)7fcS£(Mi)^ 

xe(*-* 2 )e(*-*i)e(t 2 -ti)] , 

where, as noted before k\ = k + q\ + q 2 and A; 2 = k + q 2 . 

We now expand the above traces in powers of q\ and q 2 which corresponds to performing 
a spatial derivative expansion. Furthermore since we are confining ourselves to interactions 
with static gauge fields it suffices to retain only those terms in the traces which yield couplings 
of the type 7r°e 1 - 7 diAo djA k - After considerable but straightforward algebra, the details of 
which are not particularly illuminating we find that the local effective vertex for the neutral 
pion coupling to static, classical E and B fields is given by 



S:™ eq = -Aj d A x g^(t) t ljk 7r 0+ (f,t) cUo(£) djAkix) 



(3.29) 



where the amplitude or strength of the time-dependent anomalous coupling is: 



ge 2 N ( 



6 Jt 



POO p 

dti J dt 2 J 



d 3 k 



x 



X 



1 + 



2e_d_ 

3 dk 2 



[Y k (t 2 , t) Y k (t, t ± ) + Z k (t 2 , t) Z k (t, h)} + 



2k 2 



Y k *(h,t 2 ) 



Y k (t, tl ) A x -(t 2 ,t) + r;(t,t 2 ) JLx^( tl ,t)+ 



-Zk(t,h) -^X£(t 2 ,t) - Zl(t,t 2 ) -^X^*(t 1} t) 



dk 2 



G(t - h)Q(t - t 2 )+ 



(3.30) 
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+ {X k {t 1 ,t 2 ) 1 + 



2k^_d_ 

3 dk 2 



-Y fc (t 2 ,t)n*(Mi) + z k {t 2 ,t) zi{t,t x )\ 



2k 2 



Y^tM) -^X k (t 2 ,t) + Y k *(t,t 2 ) A X -(t 1; t)- 



+Zt(t,t l )—X+(t 2 ,t) + Z* k (t,t 2 )—X+(t l ,t) f,x^ • 

e(t-t 1 )6(t-t 2 )e(t 2 -t 1 )]. 

The variables X^, Kj., are related to the mode functions f k via the following defini- 
tions: 



X k (t,t') = f k (t)f* k (t') \S k (t)S* k (t') - k 2 ] (1 - n/ ) - f* k (t)fk(t') \smSk(t') - k 2 ] n f , (3.31) 



X+(t,t') = f k (t)f* k (t') \S k (t)S* k (t') + k 2 } (1 - n f ) + r k (t)f k (t') \s* k (t)S k (t') + k 2 } n h (3.32) 



Y k (t, = -f k (t)f* k (t') [S k (t) + S* k (t')} (1 - n/ ) - f* k (t)f k (t>) [S* k (t) + S k (t')} n f , 



(3.33) 



Z k (t, t') = f k (t)f* k (t') [St(t') - S k (t)} (1 - n f ) + r k (t)fk(t') [S k (t') - S* k (t)} n f , 



(3.34) 



Sk(t)=i^ + ga (t). 



(3.35) 



While these expressions look rather unwieldy, it is clear that the 7r°77 vertex receives 
important nonequilibrium corrections through the evolution of the condensate. In the self- 
consistent mean-field treatment advocated here, we see that the effective vertex out of equi- 
librium receives contributions from all the field fluctuations through the back-reaction. 

These are memory effects since the anomalous coupling at time t depends on the fields 
at all times t' < t. Analogous effects are known to appear in related phenomena as in the 
nonequilibrium dynamics of the condensate PJ7|,P,p7|, p^ , |30| l where pion spinodal fluctuations 
and the nonlinear backreaction are very important. Hence, just as in the previous study 
of the effective vertex in equilibrium we conclude that, unlike the anomaly equation (see 
Appendix A), the effective pseudoscalar-photon vertex receives important corrections from 
long-wavelength fluctuations and is not determined by the anomaly equation. 

We are now in position to obtain the effective equations of motion for the neutral pion 
field including the self-consistent mean field evolution as well as the effective low energy 
interaction vertex with the classical electromagnetic fields: 
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d^TT°{x, t) + M 2 (t) 7r°(f , t) = -g^yit) E{x, t) ■ B(x, t) (3.36) 
M 2 (t) = -AX vl + 4A al(t) + 4A (tt 2 (x)) (3.37) 



where the last term in (|3.37|) can be read from Eq. 



This equation displays several important physical processes: i) early time spinodal in- 
stabilities resulting from the negative contribution to the mass term, just as in the more 
conventional treatments of formation of DCC []6|,[7|,P J27|J29[| , ii) an inhomogeneity induced by 



the classical quasi-static electromagnetic fields as a result of the anomalous coupling but 
with a vertex that is a strong and rapid function of time. This vertex must be obtained 
self-consistently from the coupled set of mean-field equations (|3.6| )- (|3.16| ) by inserting the 
mode functions in Eq. ( |3.31| ). Although we do not intend to offer a numerical analysis of 
the resulting equations in this article, the nonequilibrium equation of motion ( |3.36|) has the 
potential for providing an enhancement of electromagnetic isospin breaking effects through 
the combination of spinodal instabilities and nonequilibrium time dependence of the pseu- 
doscalar vertex. This enhancement could lead to an experimentally observable signal in 
the neutral to charged pion ratio. We postpone to a forthcoming article the full numerical 
study of these equations and an assessment of the potential phenomenological impact of the 
nonequilibrium dynamics. 

IV. DISCUSSIONS, CONCLUSIONS AND OUTLOOK 

The axial anomaly displays a wealth of fascinating phenomena at finite temperature 
as well as in out of equilibrium situations, perhaps even more than at zero temperature. 
What we have done in this work is to use the techniques of real time quantum field theory 
to understand the relationship between the axial anomaly and the 7r° — > 77 amplitude 
and develop a a systematic description of the in-medium equilibrium and nonequilibrium 
corrections to the pseudoscalar-photon vertex. 

The real time equilibrium calculation allowed us to make contact with some previous 
results. In particular, we obtained agreement with the results in |T7j by finding that the 
usual perturbative calculation of the triangle diagram in real time and finite temperature in 
terms of bare quark propagators is afflicted by infrared divergences in the strict chiral limit 
of vanishing quark masses, and pinch singularities in the local limit of the vertex. The latter 
singularities for vanishing external four momentum are a manifestation of the propagation of 
undamped on-shell intermediate states in the medium. Both types of divergences arise from 
the soft momentum region for the internal quark lines. In this region the quark propagators 
must be resummed including the HTL contribution from scalars and gauge bosons. The pinch 
singularities, however, are not cured by the HTL resummation but require going beyond 
the HTL limit and accounting for the on-shell width of the collective excitations. These 
conclusions are still somewhat tentative, since we have not included the vertex corrections 
whose contributions can be significant (see for e.g. |}47||). 
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After the infrared and pinch singularities had been smoothed out by HTL resummation 
and the inclusion of the quasiparticle width, we found that in the strict chiral limit, the 
anomalous vertex vanishes when m q — > in agreement with the original results of Pisarski 
13| , [T4]j3T| and the suggestion by Gelis [[17J . 

We next focused on the description of nonequilibrium dynamics during a quenched chiral 
phase transition described by a constituent quark model with a scalar-pseudoscalar sector 
described by the linear sigma model. The dynamics of a quenched chiral phase transition 
is studied in a self-consistent mean-field theory as advocated in many treatments of DCC 
formation H-fl'lH but including the dynamics of the quarks. The nonequilibrium evolution 
of the chiral condensate introduces a time dependent mass for the quarks. We obtained 
an analytic expression for the quark propagators in the mean field limit and used these to 
compute the triangle diagram out of equilibrium. The effective local pseudoscalar-photon 
vertex out of equilibrium is then extracted for quasistatic, classical electromagnetic field 
configurations as argued to be relevant in peripheral ultrarelativistic heavy ion collisions. 
An interesting aspect of the nonequilibrium calculation is that its treatment as an initial 
value problem has eliminated the infrared and pinch singularities found in the equilibrium 
situation. 

From the nonequilibrium vertex we constructed the nonequilibrium equations of motion 
for the neutral pion in the mean field approximation, including the effective pseudoscalar- 
photon vertex. This equation of motion revealed potential enhancements of electromagnetic 
isospin breaking effects through spinodal instabilities and the nonequilibrium anomalous 
vertex. Thus the triangle diagram and the effective non-equilibrium pseudoscalar-photon 
vertex is sensitive to long-wavelength fluctuations of scalars and pseudoscalars through back- 
reaction effects. 

One of the important observations (originally made in fl4jl ) of this article is that whereas 
the anomaly equation and the axial Ward identity are independent of the model and the 
low energy sector of the theory, the pseudoscalar-photon vertex in the medium in or out 
of equilibrium is non-universal, model dependent and very sensitive to low energy physics. 
The next calculation to do would be a numerical study of the effective equation of motion 
for the neutral pion, with an eye towards answering the question addressed in [p0[| : can the 
anomalous vertex enhance the formation of neutral pion DCC regions. This would have 
fascinating phenomenological consequences. 
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APPENDIX A: THE ANOMALOUS WARD IDENTITY OUT OF EQUILIBRIUM: 



The fact that the axial Ward identity is unmodified by thermal corrections in a heat 
bath at finite temperature was explicitly shown in [[TJJ. In this appendix we explicitly show 
via perturbation theory at 1-loop that the axial Ward identity and the anomaly equation 
are not modified in or out of equilibrium. This is a non-trivial statement in the real time 
formulation since there are, as we will see below, different contributions to the triangle 
diagram and anomaly equation arising from the different branches in the CTP contour of 
the generating functional. 

We will use the point-splitting approach to derive the anomalous conservation law but 
implemented in real time. 

The classically conserved axial current in the LSM is 

& = VWT°V + - n a d,a , (Al) 

where the superscript (a) denotes the isospin component under consideration. The quantity 
of interest is the a = 3 component of the current which exhibits the anomalous behaviour. In 
the nonequilibrium CTP formulation all the operators reside on the CTP contour, carrying 
+ and — labels depending on which part of the time contour they belong to and we will 
compute within perturbation theory the expectation value of 

WST>- (A2) 

which is defined on the forward contour. In the subsequent steps we will omit the isospin 
index since it is clear which component we are focussing on. The expectation value above is 
of course a priori ill-defined because of the usual short- distance singularities in the operator 
product expansion. We therefore regulate this composite operator by point-splitting in a 
gauge invariant manner as follows, 

(Pj&\x))=^{Pj£\x,e)) } (A3) 
where we have defined a new operator 

(A4) 

i> + {x + e) llll5 r 3 e ieQ ^ +CA{zydz ^ + (x) + a + (x)dy + (x) - n 0+ (x)d II a + (x) . 



The Wilson-line has been introduced to render the point-split operator gauge- invariant. This 
is the analogous to the treatment of the axial anomaly in spinor QED for the vacuum 



The gauge-fields can be treated simply as background fields without loss of generality. 
Using the equations of motion it is now easy to see that 

(d^(x, e)> = h(7r 0+ (x)) + te(^ + (x + e)Yl,Qr^ + {x))F^{l + 0(e)) (A5) 
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The correlator on the right hand side of Eq. ( |A5| ) has short distance singularities as e M — ► 0, 
which we now extract via a one-loop calculation. The question that we want to answer is 
whether the nonequilibrium path integral introduces any new divergent contributions to the 
Ward identity. At one-loop we find 



<V> + (x + t)^Qr^ + {x)) = -ie^- J d 4 z {Tr [S ++ (z, x + e)^ 5 S ++ (x, z)<y 



+ 



(A6) 



-Tr 



■(z, x + e)n 5 S + -(x, z)>f] } A v {z) 



Rewriting Eq. 
we obtain 



in terms of the spatial Fourier transforms of the Green's functions 
d 3 k 



6 J "' J (2tt) 3 'J (2tt) 3 

{Tr[5+ + (t',t + eV75^ + +(t,t')7 



(A7) 



Tr 



S>(t',t + e°)^ l5 S<(t,t') k+q Y}} 



There are two distinct contributions to the Ward identity out-of-equilibrium. The first 
one is obtained by contractions with fields on the forward time contour and involves the 
usual time-ordered propagators. The second term which arises from contractions with fields 
on the backward time contour is the new contribution that appears in the nonequilibrium 
case. Non-trivial modifications of the Ward identity can only arise from terms that diverge 
as e — > 0. But the latter are short distance divergences and therefore we need only to 
consider the UV behaviour of the integrand in Eq. (A~7). As \k\ —>■ 00, go~ 2 (t) ga(t) and 
the temperature % can be ignored and the propagators in Eqs. ( |3.11[) , (|3.12|) and ( [3.14|) 
take on their usual zero temperature short distance forms, which is simply dictated by the 
singularities of free-field theory in the operator product expansion of fields. Hence the UV 
portion of the first term can be rewritten in terms of the Fourier transformed fields as 



6 J (2vr) 4 



(2vr) 4 k 2 (k + q) 2 



(A* 



which, after evaluating the divergent integral, yields 



-ieN c f e a 
12tt 5 



^f) d p A v {x)e a ^. 



(A9) 



So the contribution of this term to the Ward identity is then 
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lim 

c— 



e„e^\ e 2 N r 



24vr 2 



Fa F e a ^ v 



e 2 N c 

967T 2 



F^yFape 



liva/3 



(A10) 
(All) 



We now isolate the ultraviolet (the leading terms in the e — > limit) contributions to the 
second term in Eq. ( |A7|) which is due to the presence of the backward time contour 



6 

x 



ip N 

!ffl£( 4l - e ^) x (A12) 



j4 „ J3t, ik-e it(\k+q\+k) f 



(2vr)4 ^' > J (2tt)3 Ak\k + q\ 



6 



(2ir)" ' J (27r)3 4 fc|fc + <3l 

Clearly the fc- integral is cut off at large momenta by the delta function and is finite in the 
limit e — > 0. Because of the presence of e u multiplying the expression (see Eq. flA"5f) ), this 
term does not survive the e — > limit, hence there will be no further anomalous contributions 
to the Ward identity other than Eq. ( |A11|) . Therefore the anomalous Ward identity is in 
fact unchanged by equilibrium finite temperature or the nonequilibrium state, as expected: 



p 2 N 

WFW) = H^ 0+ ) - ^ F ^ vaP - ( A13 ) 



That there are no corrections in or out of equilibrium should come as no surprise since 
the anomalous contributions to the axial Ward identity are generated by the ultraviolet 
(short distance) properties of correlators and should not be affected by either temperature, 
finite density or the time dependence of background fields or long-wavelength fluctuations. 
However, since a detailed calculation in a medium in real time in or out of equilibrium 
was missing, this analysis provides a reassuring generalization of the vacuum result. The 
calculation of the expectation value for the backward branch follows the same steps and 
leads to the same results. 
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APPENDIX B: NONEQUILIBRIUM FERMION GREEN FUNCTIONS 



Below we outline the construction of the nonequilibrium Green's functions for the 
fermions in the presence of the time-dependent background condensate cxo(t) ||28|| . We con- 
sider the solutions of the time dependent Dirac equation 



[i$ - ga (t)]ij)(£,t) 
Writing the four independent solutions as 







(Bl) 



V (1 ' 2) (x,t) 

the mode functions are found to satisfy 



d -> r i \ 

d _ r> . . 

z 7o^ + 7 • k - ga {t) 



ik-xj r(lj2) 



ur>(t) 



-ik-x-r a(1i2) 



o. 

. 



(B2) 
(B3) 



It turns out that it is convenient to write the spinors as 



U^\t) 



v^\t) 



d _ r> . . 

^o-j: - 1 ■ k + gcr (t) 
dt 

d - r / \ 
no-jr + 1 ■ k + ga {t) 
dt 



fk(t)u^ 



(B4) 
(B5) 



with 

u (l,2) j v (l,2) 

the spinor eigenstates of 70 with eigenvalues +1, —1 respectively. The 
functions fk{t), gkif) obey the second order equations 



dt 2 
' £_ 
d^ 



+ k 2 + g 2 a 2 (t)-iga (t) 
+ k 2 + g 2 a 2 {t) + iga {t) 



fk{t) 
9k{t) 



■ 

: 



(B6) 
(B7) 



We now need to append initial conditions. We will consider the situation in which the system 
was in equilibrium at a temperature Tj for all times t < tj with the expectation value of the 
scalar field being cx(0) and d(0) = 0. Thus the fermion mass is constant and given by go"o(0). 
We can now impose the condition that the modes fk(t), guit) describe positive and negative 
frequency solutions for t <ti and, normalizing the spinor solutions to the Dirac equation to 
unity, we impose the following initial conditions that describe particle and antiparticle states 
for t < tj < 
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f k (t < 0) = = , (B8) 

'2uj k (uj k + ga (0)) 



e iu k t 



g k (t < 0) = , (B9) 

^2u k (2u k + ga {0)) 



u k = Vfc 2 + pV 2 (0). (BIO) 
Equations (|B7|) with these boundary conditions imply that 

g k (t) = f k (t). (Bll) 



1. Initial temperature Tj = 0: 



The necessary ingredients for the zero temperature fermionic Green's functions are the 
following 



a=l,2 



o=l,2 



In the standard Dirac representation for the 7 matrices, we find 
Sok(ti t ) = 

-ifk(t)fW) [Wk(t)y -y-k + ga Q (t) 

-if k (t)fk{t') [Wmio -J-k- ga {t) 
./*(*) 



l + 7o 
2 

1 -7o 



W* k (t')ya-y-k + g*o(t') 
m(t'ho-i-k-ga (t') 



m(t) 



f k (t) 



(B12) 
(B13) 

(B14) 
(B15) 

(B16) 



With these basic results in place, the nonequilibrium fermionic Green's functions can 
be constructed using Eqs. ( |3.7| ), fl3.7|) , ( |3.9| ) and Q3.10]) . The fact that the equality 
Tr So~ k + (t, t) = Tr S ok T(t,t) = Tr 5^.(£, t) = Tr S^ k (t, t) is satisfied (where the trace is over 
Dirac indices) provides an important check. The mode functions f k (t) satisfy yet another im- 
portant property which is a consequence of probability conservation and can also be checked 
explicitly pi, 



iti(t)+g<T (t)ft(t) if k (t)+ga (t)f k (t) + k l ft(t)f k (t) 



(B17) 
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2. Initial temperature Tj 7^ 0: 



When the initial temperature is not zero, the CTP contour extends into the complex 
plane with a leg that runs from t to t — if3.. The fermionic fields and Green functions then 
satisfy antiperiodic boundary conditions: 



S>(t -zf3,t) = -S<(t ,t). 



(B18) 



Incorporating these boundary conditions then leads to the following formulae for the nonequi- 
librium fermion propagators at finite initial temperature: 



S?(t,t') = S>M,t') [1 - n f M\ - S<M,t')n f (u; k ) , 



(B19) 



^i^') = S&M [1 - n/K)] - S>^t : t')n f (u; k ) , 



Ok y 



(B20) 



k 2 + g 2 a 2 (0) , 



n f (u k ) 



e Poj k _|_ I ■ 



In equilibrium, i.e. when cr is a constant function of time, the fermion propagators are: 
' r e -*"*(*-0( # + rrty)(l - n/K)) + e^-^o^ - m^To^/M 



2u; fc 



2oa 



(B21) 
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